CERN-TH. 6323/91 
November 1991 



On 
Os 



On 



> 

o 

(N 



OS 



X 



Universality and Non-Perturbative Definitions 
of 2D Quantum Gravity from Matrix Models 



J.Luis Miramontes 



o 

Q \ Theory Division, CERN 



CH-1211 Geneve 23, Switzerland 



and 



Joaquin Sanchez Guillen 



Department of Theoretical Physics 



^| University of Lund 



Solvegatan lJ^A 
S-22362 Lund, Sweden 

ABSTRACT 

The universality of the non-perturbative definition of Hermitian one-matrix 
models following the quantum, stochastic, or d = 1-like stabilization is discussed in 
comparison with other procedures. We also present another alternative definition, 
which illustrates the need of new physical input for d = matrix models to make 
contact with 2D quantum gravity at the non-perturbative level. 
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1. INTRODUCTION 



The double-scaling limit of the matrix models provides a regularized version of the sum 
over arbitrary surfaces and, consequently, of 2D quantum gravity 1,2 . In the case of the 
Hermitian one- matrix models the result, by now well-known, is that there is an infinite 
number of local operators, Ok, k = 0,1,..., which realize a KdV flow structure 2 ' 3 , and 
whose simplest operator is the puncture Oo = P. Moreover, the specific heat of the model 
satisfies a differential equation known as the "string equation", whose scaling solutions in 
the cosmological constant are the multicritical models 4 . Introducing the couplings of the 
operators Ok, t k , the string equation is 

oo 

^2(2k + l)t k R k [f\ = -T, (1.1) 
fe=i 

where T = to is the renormalized cosmological constant, / = <i 2 lnZ/<iT 2 = (PP) is the 
specific heat, R k are the Gelfand-Dikii polynomials that define the KdV flows 

J£ = (PPO k ) = R' k (1.2) 

and primes denote derivatives with respect to T. The k tli multicritical model, which is defined 
by the asymptotic behaviour / pa T x l k + . . . for large positive values of T, corresponds to 
tj = (2k + l) -1 ^. The simplest one realized in the matrix model is pure gravity (k — 2), 
whose string equation is the Painleve I equation 

f ~ \f" = T. (1.3) 

In any of these models, the contribution of the surfaces with an arbitrary fixed genus to 
the specific heat is obtained by expanding the solution of the string equation in powers of the 

2fc + 1 

string coupling constant, which is related to the cosmological constant, g s tr oc T ~ . It was 
also expected that the string equation, which, strictly speaking, is only valid order by order in 
the genus expansion, would define 2D quantum gravity at the non-perturbative level as well. 
This would require to fix the string equation's boundary conditions. Nevertheless, although 
it is possible for the odd-A; models 5,6 , the even-k ones present additional complications. In 
the particular case of pure gravity, it has been argued that no real solution of the Painleve I 
equation can be compatible with the loop equations (or Schwinger-Dyson equations) of the 
matrix model 7 . Furthermore, the behaviour of the (KdV) flow between the (well defined) 
k = 3 and k — 2 models indicates the non-perturbative instability of pure gravity 8 . 
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The definition of the models with even-A; involves matrix models whose potentials are 
not bounded from below, and this seems to be the origin of the problems 5 . Consequently, the 
two standard procedures to stabilize field theories with unbounded actions have been tried 
in this context. The first method is the analytic continuation of the dominant term of the 
potential from positive to negative sign, while a simultaneous deformation of the contour of 
integration on the matrix eigenvalues is performed 9 . Unfortunately, this method leads to a 
non-perturbative imaginary part in the specific heat, which makes the result unphysical 7 ' 10 ' 11 . 
Using the second generic method, inspired in stochastic quantization and proposed in general 
by Greensite and Halpern 12 , the ill-defined (d = 0) matrix model is formulated like a well- 
defined (d — 1) quantum mechanical system 13 ' 14 ' 15 . Following this method, the specific heat 
remains real, in contrast with the result of the analytic continuation. 

In addition to its relevance in the context of the matrix models, the study of the 
stochastic-like or quantum mechanical method is interesting by itself, as a non-trivial exam- 
ple of the general stabilization mechanism proposed in ref.[12]. So far, the universality of this 
regularization has not been fully understood. In fact, all the previous papers concentrate 
on the definition of pure gravity using the simplest potentials of degree 3 or 4, either in the 
WKB approximation 13 ' 14 , or in numerical computations 15 ' 17 . The main reason for this lack 
of generality is that almost all the calculations rely on the fact that the quantum mechanical 
system is an ideal Fermi gas for the simplest potentials of degree less than 4. In the next 
section 2, we show that it is possible to recover the ideal Fermi gas picture for an arbitrary 
potential through the use of a mean field approximation (a la Hartree-Fock), which agrees 
with the semiclassical limit of the matrix model given by the loop equations. This result 
ensures the universality of this regularization procedure for any given multicritical model 
(not only pure gravity), and provides a connection between the quantum (Fokker-Planck) 
potential and the semiclassical density of eigenvalues. Using this connection, we show in sec- 
tion 3 that the dominant non-perturbative effects in the quantum mechanical regularization 
of pure gravity are given by a metastable instanton (bounce) whose lifetime is universal and 
of the size of the non-perturbative ambiguities in the specific heat. These results ensure, in 
particular, that the method of Greensite and Halpern is compatible with the double-scaling 
limit. 

a In this paper, we shall not be interested in the supersymmetric aspects of the original proposal of 
Marinari and Parisi 13 ' 16 . 
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The problem of pure gravity and of all the even-A; models is that their genus expan- 
sions for the specific heat are not Borel-summable; hence, they are not well defined at the 
non-perturbative level 2,18 . Therefore, one should not be surprised to have many different 
"definitions" of pure gravity if the agreement with its genus expansion is the only select- 
ing criterion. In fact, a third consistent definition has been proposed recently by imposing 
the non-perturbative agreement with the KdV flow structure of eqs. (1.1) and (1.2) 19 . 
This proposal provides a real specific heat too, which is different from that defined by the 
stochastic-like regularization 20 . Obviously, some new physical insight into non-perturbative 
2D quantum gravity is required to select which is the right definition. A good example of a 
true physical constraint is the unitarity of the ^-matrix in the d — 1 matrix model, where 
a clear spacetime interpretation is available, in contrast with the d = case 21 . In fact, 
in d = 0, the agreement with the KdV flow structure could be a good criterion, but the 
matrix models provide the KdV structure only at the perturbative level (genus by genus), 
and it can be broken through non-perturbative contributions. To conclude the paper, we 
discuss in section 4 the solutions of the Hermitian and anti-Hermitian 1-matrix models with 
generic potentials that illustrate this last point. Furthermore, we argue that the solution 
of the anti-Hermitian model is the simplest way to accommodate the complex solutions of 
the Painleve equation that have been considered in the analytic continuation method, and 
provides another consistent definition of pure gravity. This shows, once more, that the real 
problem is our ignorance about quantum gravity at the non-perturbative level, and that the 
matrix models alone do not solve these ambiguities. 

2. WKB APPROXIMATION AND THE 
SEMICLASSICAL MATRIX MODEL 

In the stochastic-like regularization method, it has been proved that the perturbative 
expansion of the equivalent quantum mechanical system is the same as that of the original 
(unstable) model provided that they agree at the first order. In the case of pure gravity, 
this has already been checked using the lowest-order potentials, where the quantum system 
is an ideal Fermi gas of N particles. In this section, we show that the quantum mechanical 
regularization reproduces, in the WKB approximation, the semiclassical limit of the matrix 
model not only in the case of pure gravity, but for any multicritical model. The main 
difficulty is the fact that the potentials with degree > 4 induce interactions between the 
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fermions, and the Fermi gas is no longer ideal. Nevertheless, a mean field approximation can 
be used to decouple the Fermi gas in the WKB approximation. 

Let us consider the d — Hermitian one-matrix model, whose partition function is 

Z= f D(f>e-P v ^\ (2.1) 



where is an N x N Hermitian matrix, and V(<f>) is a generic potential of degree L 

L 

V(0) = 5>"Tr(0 n ). (2.2) 

Following refs. [14,15], this model can be formulated like the ground-state of a d — 1 quan- 
tum mechanical system, whose Hamiltonian is the (positive semi-definite) Fokker-Planck 
Hamiltonian 

a o2 ay ay o p[2y 

H FP = Tr(P 2 ) + W FP with Py = -i— and W FP = '^r^-£- ~ g ' , , , • (2.3) 

d(f)ji 4 d^ij dcpji 2 d^ijd^ji 

This means that the expectation value of a generic operator Q{4>) in the matrix model 
corresponds to the vacuum-expectation- value (VEV) of the quantum operator Q(<f>): 

(Q) = (0| Q($) |0) = y P0^(0)g(0), P FP ^o(0) = ^o^o(0), (2.4) 

where Eq and x I / o(0) are the energy and the wave function of the ground-state. Eq = 
corresponds to the case of potentials bounded from below, whose matrix models are well 
defined, and ^o(0) — exp(— j3V((fi)/2). Otherwise, if Eq > 0, eq.(2.4) defines (Q) in terms 
of the true ground-state. Notice that the extra dimension introduced is just an auxiliary 
degree of freedom without interpretation within the matrix model. In fact, the correlators 
of operators taken at different times in the quantum theory do not have any meaning in 
terms of the matrix model. At this point, it is worthwhile to mention that the ground-state 
energy Eq is not the partition function of the d = matrix model. In fact, this method only 
provides closed expressions for the correlators. Nevertheless, it is easy to obtain directly the 
derivatives of the matrix model partition function with respect to the cosmological constant, 
i.e., the correlation functions of the puncture operator P. They are the universal scaling 
parts of the correlation functions of Tr0 n , for any finite value of n (see eq.(2.20)). Therefore, 
one could add a perturbation to the Fokker-Planck Hamiltonian, Hpp — > Hpp + JTr(0 n ), 
such that the derivative of the perturbed ground-state energy Eq(J) with respect to the 
source J, evaluated at J = 0, gives the correlation function of the puncture 17 . 
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Following the well-known techniques of ref . [22] , it is natural to make a change of variables 
to the eigenvalues {Aj} of the matrix 0, by introducing the effective ground-state wave 
function 

*S // ({A i })=II( A *- A i)*o(0), (2-5) 

i<j 

which is totally antisymmetric, and describes a gas of N Fermi particles. In general, the 
Fokker-Planck potential, Wpp, can be splitted into its diagonal (D) and non-diagonal (ND) 
parts in terms of the eigenvalues 

w FP = wj$ + 



(2.6) 



< D) = -f EIE^E^ 

\n=4 s=0 



iV / L n-4 



n— 3— s 



where X = N//3 = e 70 is related to the (bare) 2D cosmological constant in the usual way. 
Obviously the ND piece does not vanish if L > 4, and the Fermi gas is not decoupled. 
Nevertheless, in the semiclassical WKB limit ((3 « h^ 1 — > oo), a mean field approximation 
(a la Hartree-Fock) may be performed to decouple the system. We show below that 

Tr(0 fc )Tr((^) « JV (w fc Tr(0 p ) + cu p Tr(0 fc ) - iVw fc w p ) + • • • (2.7) 

where the normalization is fixed by 

<Tr(0 fc )Tr(0^)) c » iV 2 (^ + 0(1/ N)) and co fc = ^<Tr(0 fc )) c , (2.8) 

is consistent with the semiclassical limit of the matrix model. Under this approximation, the 
quantum mechanical system becomes an ideal Fermi gas of N particles whose Hamiltonian 
is 

N o2 o2 

H FP ^J2 hFP ( x i^ with hFP ^ = ~W + T f/Fp(A) ' and 

i=l 

(L-2 I L A L t-2 

E E wwj-i-2 a 1 + (? 2 - 2 E ^ E ^-i^i-j-i 

(2.9) 

The effective potential ?7t?p is bounded from below, and the one-fermion Hamiltonian has a 
well-defined discrete spectrum, h FP ip n (\) = 4-e n ?/> n (A). The ground-state wave function is 
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the Slater determinant of the first N eigenfunctions, and the vacuum expectation value of a 
generic operator Q(4>) is 



(O|Q(0)|O) = J [dA]det z (^_i(Aj))g({A}), i,j = l,...,N, (2.10) 
to be compared with the corresponding expression in the matrix model 

<W)> = ! J[d\]det 2 (e-^)Aj- 1 ) Q({A}), i,j — l,...,N. (2.11) 

All the relevant information about the d = 1 quantum mechanical system is contained 
in the particle density, p(A, e) = (A| 5 (h FP — ^-ej |A), whose normalization fixes the Fermi 
energy e F : 

e F 

N = j d\ J de p(A,e). (2.12) 

— oo 

The energy integral of p(A, e) provides the quantum-mechanical version of the matrix model 
semiclassical density of eigenvalues 

u(X) = i J de p(A,e), X = j = J dX «(A), (2.13) 

— oo 

which, in the WKB approximation, is 

u WKB {\) = ^^/e F -U FP {\) 9(e F - U FP (X)) 
i <Trg(0)>= J d\ Q{\)u WKB {\). 



(2.14) 



In order to compare the WKB result, eq.(2.14), with the matrix model, we shall use the 
semiclassical limit of the Schwinger-Dyson loop equations 7 . In the semiclassical (planar) 
limit, the generating function of monomial expectation values 

satisfies the Schwinger-Dyson equation 

L-2 / L \ 

F(V? ~ V(p)F(p) + X E 39j P* = 0, (2.16) 

i=0 \j=i+2 J 

where the "constants of integration" oj^ have been already defined in eq.(2.8). The solution 



6 



of this equation is 

F(p) = JdX = \ (V(p) - y/ZQj) 

l-2 ( l \ (2.17) 

Hp) = (V(p)) 2 -4xJ2 E v\ 

i=0 \j=i+2 J 

and, for a given multicritical model, the "constants of integration" are fixed by the condition 
that the imaginary part of F(p) defines a proper semiclassical density of eigenvalues 7 ' 22 . 
Obviously, A(p) is a polynomial in p and all the singularities of F(p) will be the branch 
cuts of y/A(p). Therefore, A(p) has to satisfy the following constraints 7 ' 22 ' 23 : (i) A(p) must 
have only real zeros in the complex p-plane, and (ii) A(p) cannot have three consecutive odd 
degree zeros. Under these conditions, the semiclassical density of eigenvalues is 

u sc (X) = -Im F(X) = ^y^AjXj 6>(-A(A)). (2.18) 

7T Z7T 

The branch cuts of F(X), i.e., the intervals between odd degree (real) zeros of A(A), are the 
"bands" on which u sc (X) has support. 6 

Therefore, under the above mentioned restrictions, the comparison between eqs.(2.17), 
(2.18) and eqs.(2.9), (2.14) shows that the WKB limit of the d = 1 Fokker-Planck Hamil- 
tonian, with the mean field approximation of eq.(2.7), agrees with the semiclassical limit of 
the matrix model if the Fermi energy is 



e F 



(L) LI i-2 



! x =92 + 3^1 + E ( Ui ~ 2 + \ Yl (^j-i^i-j-i) m- (2-19) 
i=4 \ j=2 J 

Notice, eq.(2.8), that the "constants of integration" are, in fact, correlation functions. 
Accordingly, the universal piece of the Fermi energy in the double-scaling limit can be ex- 
pressed in terms of the expectation value of the puncture operator. In the case of even 
potentials, the precise relation between Uk and the puncture operator is, for the k th model 2 , 

Nu 2s = (Tr(0 2s )) w n.u. + s^)/r^(P), (2.20) 

where n.u. stands for the non-universal part. Therefore, the dominant universal piece of the 



b All the multicritical models of refs.[l,2] correspond to single-band configurations, which means that 
A(p) has exactly two single (odd degree, in general) real zeros in the complex p-plane. 
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Fermi energy is 

. + («f (< iy! ! 2)! (^))r^<») + ■■■ (2.2D 

In fact, this is the most direct way to obtain the critical exponents of the multicritical models 
within the quantum mechanical formalism. Moreover, the scaling factor when f3 — > oo will 
dictate the relevant piece of the Fokker-Planck potential in the double-scaling limit. 

Our result agrees with, and generalizes, previous results obtained for pure gravity with 
the simplest potentials of degree three 13 ' 14 ' 15 ' 17 ' 24 and four (even) 25 , showing that the critical 
behaviour of the Fokker-Planck Hamiltonian is precisely that of the matrix model. Therefore, 
in the semiclassical approximation, it is possible to describe any Hermitian one-matrix model 
as an ideal Fermi gas of N particles, whose potential is given by the semiclassical density of 
eigenvalues 

U FP (X) =e F - [2nu sc (X)} 2 . (2.22) 

Notice that this relationship formally holds only when the above mentioned constraints 
on the zeros of A = Upp — e F are satisfied and u sc is well defined. Nevertheless, the 
quantum mechanical system, and u WKB , is defined even when this is not the case. Such 
quantum mechanical configurations arise when the naive ground-state wave function, \l>o = 
exp(— PV/2), is not normalizable in the WKB approximation either. Therefore, they do 
not correspond to any semiclassical configuration of the matrix model, and they contain 
non-perturbative information about the stabilization mechanism. In fact, their study has 
shown that the quantum mechanical formulation is not compatible with the KdV flows at 
the non-perturbative level 20 . 

3. UNIVERSAL NON-PERTURBATIVE BEHAVIOUR 

In the previous section, we have shown that the WKB limit of the quantum mechanical 
formulation agrees with the semiclassical limit of the matrix model for an arbitrary potential. 
If the stabilization mechanism is compatible with the double-scaling limit, this result ensures 
that the quantum mechanical formalism preserves the universality of the multicritical mod- 
els. In this section, we show that the non-perturbative contributions are indeed universal 
in the case of pure gravity, and that they agree qualitatively with the ambiguities expected 
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in the specific heat. The size of these ambiguities is obtained through the linearization of 
the Painleve equation, eq.(1.3). The result is that the difference between any two func- 
tions whose genus expansion is the pure gravity genus expansion should be propotional to 
exp ^— i^lT 5 / 4 j when T — > oo 18 ' 26 . We show below that the typical non-perturbative 
exponential suppression corresponds to a universal mestastable instanton (bounce) in the 
quantum mechanical picture, which is a natural relationship in quantum mechanics when 
the perturbative expansion is not Borel-summable. 

To identify the universal features of the Fokker-Planck potential, the main tool we shall 
use is the connection with the semiclassical density of eigenvalues, eq.(2.22), which allows 
the universal behaviour of the potential to be derived from that of the semiclassical density. 
From now on, we shall restrict ourselves to the case of the standard multicritical models 
defined with even potentials, but the results should be easily extended to the general case. 
Using the standard orthogonal polynomial techniques 27 , we can compute the semiclassical 
density of eigenvalues directly from the generating function of monomial expectation values, 
eqs.(2.16), (2.17), and (2.18) 2 ' 23 : 

N-l 



N-l 



n=0 
X 

1 /" j 1 

ax 



(3.1) 



iy 

^ \/P 2 -4i?(n) J y/p'-ARix)' 



where R n is the coefficient in the recurrence relation of the polynomials, XP n = P n +i + 
R n P n -i, x = n/P, and R(x) is the limit of R n when (5 — > oo. The imaginary part of eq.(3.1) 
provides the semiclassical density of eigenvalues 22 

x 

u sc (p) = -ImF(p) = - I dx . 1 6{AR{x) - p 2 ). (3.2) 

7T J JAR(x) - V 2 







Next, we make the change of variables x = W(R(x)), with W(R) defined in terms of the 
potential as 
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and the final expression for the semiclassical density is 

R(X) 



p 2 /4 P 

= ^-V e F- U FP {p)6{e F - U FP (p)). 

271 



(3.4) 



For the k th model, the double-scaling limit corresponds to (5 — > oo, with c 

W(R) = 1 - (1 - R) k , X = 1- (3-^TiT, and R = 1 - /T^/(T). (3.5) 
Therefore, the critical behaviour of use is 

v/4i?(X)-p2 

? (p) = J* / ^(4-P 2 -y 2 ) fc_1 . (3.6) 



Notice that Upp(X) has two roots of order 2/c — 1, A = ±2, at the critical point X = R = 1. 
Thus, the resulting Fokker-Planck potential has degree 4k — 2, and it corresponds to a 
matrix model whose potential has degree L = 2k, i.e., the canonical representative of the 
multicritical /c th model 2 . 

The Fokker-Planck potential of pure gravity, k — 2, is 

C/ FP (A) - e F = - (A 2 - 4R(X)) (A 2 - 6 + 2R(X)) 2 . (3.7) 
9 

This potential has an absolute minimum at A = 0, and two relative minima at A 2 = 6 — 
2R(X); obviously, the Fermi energy is just the value of the potential at these relative minima. 
The perturbative expansion around the absolute minimum provides the WKB expansion, 
which reproduces the 1/N expansion of the matrix model 12 . Nevertheless, we are interested 
in the double-scaling limit, which corresponds to /3 — > oo. We have already obtained the 
scaling behaviour of the Fermi energy in this limit, eq.(2.21), which is like /3~ 6//5 for pure 
gravity. This behaviour corresponds to the scaling of the potential around the secondary 
minima, i.e., around the Fermi energy, which is^ 

6 16 

U FP (X) w ±/T s — (4z 3 - 3Tz) +■■■ (3.8) 
9 

2 

where z = /3s (A =F 2). Notice that, in this limit, the potential becomes unbounded from 



c As is usual, we have fixed the critical values X c = R c = 1, but the result is independent of this choice 28 . 
d Even though we have obtained this expression starting from even potentials, it can also be derived 
from the cubic potential 15 , as expected from its universal character. 
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below. However, this does not destroy the stabilization mechanism. In fact, the subdominant 
terms in the j3 — > oo limit are crucial to ensure that the potential has bound-states, and, 
in particular, the normalization condition in eq.(2.13). On the other hand, the universal 
properties correspond to the behaviour of the potential around the Fermi level, and only the 
energy levels close to the N th one are relevant in this limit. Because N — > oo, these levels do 
not feel the absolute minimum of the potential, which appear to be unbounded from below. 
In fact, one can check that the behaviour of the potential around A = does depend on the 
potential one chooses to define pure gravity, and hence it is not universal. 

The dominant non-perturbative contributions are related to the classical possibility that 
the fermions at the Fermi level could be at the secondary minimum of the potential, falling 
later into the main (unbounded) well due to quantum effects (tunneling). This configuration 
corresponds to a metastable instanton, also known as a bounce 29 , which decays through 
barrier penetration effects. The lifetime of the bounce is proportional to the imaginary 
part of its energy, which can be computed from the Euclidean partition function in the 
semiclassical approximation 29 ' 30 . 6 The Euclidean partition function corresponding to the 
Hamiltonian of eq.(2.9), h FP (\) = + ^U FP (X), can be written as the L — > oo limit of 

the following functional integral 

Tr(V Li/ )= J [dq(t)}e~ S[q{t) \ (3.9) 

q{-L/2)=q{+L/2) 

where S[q(t)} is the Euclidean action 

+L/2 

S[q(t)]= J dt(^q\t) + ^U FP m))- (3-10) 

-L/2 

Tr(e~ LH ) is given by ^2 ex P ( — LE n ), where the sum extends over the whole spectrum. 
Consequently, in the L — > oo limit, the partition function is dominated by the ground- 
state energy. Nevertheless, the barrier penetration effects induce an imaginary part in the 
energy value of the bounce, Ef, c , and, then, in Tr(e~ LH ). The dominant contribution can 
be obtained by computing the imaginary part of the partition function in the semiclassical 

e Notice that this imaginary part cannot be identified with the one obtained in the analytical continuation 
method, because the ground-state energy of the quantum mechanical system is not the specific heat of 
pure gravity, which is real within the stochastic-like regularization. 
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B 2 

approximation. If E^ c = Ep + iT, where Ep = ^-ep, then, in the L — > oo limit, 

Im (Ti{e- LH )^ oc Im ^ e -HE F + %T)^ w LT (T LEf , (3.11) 

where we have taken into account the smallness of the non-perturbative imaginary part. 

In the semiclassical approximation, the dominant imaginary contribution to Ti{e~ LH ) 

is given by the bounce, which is a solution of the Euclidean equations of motion, ^q c = 

b 2 

4-U F p(q c ), that starts from the secondary minimum, q m , at Euclidean time L = — oo, is 
reflected in the classical turning point, qo, and comes back to the minimum at time L = +00. 
Therefore, the bounce action is 

+00 

%d = Jdt Qg c 2 (t) + ^U FP (q c ^ 

(3.12) 

= LE F + J dq^(3 2 {U F p{q)-e F ) = L Ep + S hc . 

The integration around q c in the Gaussian approximation provides the bounce contribution 
to the partition function. This calculation involves a Jacobian factor that becomes complex 
because of the turning point qo, and is proportional to L because of translation invariance 
in the Euclidean time. The final result is 

Im (Tr(e" Li/ )) oc Le^] = Le -LE F -S bc w Im ^-L(E F+ i e -^ ^ 

Therefore, T oc exp(— Sf, c ) is the imaginary part of the bounce energy, and the lifetime of 
the bounce is r^ c oc T. In the case of pure gravity, Sb c remains finite in the double-scaling 
limit 28 , and its value is 

s fc = i^rf. (3.i4) 



The lifetime of the bounce is exponentially suppressed for large values of T, and its value 
is expected to give the characteristic size of non-perturbative effects. Notice that exp(— S^c) 
is, in fact, the size of the ambiguities between functions whose genus expansion is that of 
the matrix model. Therefore, we conclude that the dominant non-perturbative effects of 
the quantum mechanics regularization of the matrix model are universal, and consistent 
with those expected for pure gravity. In fact, this checks, in a non-trivial way, that the 
stabilization procedure is consistent with the double-scaling limit because it ensures that the 
genus expansion coming from this method agrees with that of the original matrix model 30 . 



12 



For completeness, let us consider eq.(3.6) for the first well-defined multicritical model, 
k = 3, 

u sc (p) oc (AR(X) -p 2 Y {{p 2 + R(X) - 5) 2 + 5 (1 - i?(X)) 2 ) . (3.15) 

In this case, Upp(X) has only one absolute minimum at A = 0. Therefore, the non- 
perturbative contributions to the k — 3 model are not related to metastable states. This 
agrees with the expected relationship between metastability and non-Borel summability 
of the perturbative series in quantum mechanics. Therefore, we expect that there will 
be bounces in all the k-even models, and not in the fc-odd ones. This conclusion gives 
support to the identification of the instabilities in the matrix model with instanton-like 
contributions 10 ' 18 . In fact, it shows that the criterion for having or not instabilities is given 
by the Fokker-Planck potential Upp in the double-scaling limit, and not directly by the 
potential of the matrix model. 

To conclude this section, let us estimate the dominant non-perturbative contributions 
to the specific heat within the quantum mechanical formalism. The matrix model partition 
function is expressed in terms of the density of eigenvalues as 22 

f mm = p2 J dXu ( X ) (^J dnu(ii) ln\X -fi\- V(X)^j . (3.16) 

The dominant non-perturbative contributions to F MM will correspond to the dominant 
non-perturbative modifications of the density of eigenvalues, Su(X), which are related to 
the bounce and are dominant outside the classical (WKB) range. Therefore, the dominant 
modification of the matrix model is 18 



6F mm = (3 2 J dX6u ^ ^ 2 J dnu WKB {n)ln\X- n\ - V(X) 

A 

= -f3 2 J dX5u{X) J dq^U FP {q) - e F 9(q 



(3.17) 

go), 



qo 

where eq.(2.17) has been used, and go is the classical turning point at the Fermi energy, i.e., 
the limit of the classically allowed range of motion, eq.(3.4). 

The dominant contribution to Su(X) is related to the classical possibility that the fermions 
at the Fermi energy can be sitting at the secondary minimum. Therefore, one should expect 
8u(X) oc J\fS(X — q m ), where the normalization is expected to be proportional to /9 _1 , because 
it involves only the fermions at the Fermi energy 18 , and to exp(— S&c), because it is the 
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expected suppression of the non perturbative effects. In fact, we can check this normalization 
with the results of ref.[24], where the non-perturbative effects in the bound-state energies 
have been computed, and they result to be proportional to /3 -6 / 5 exp(— S& c ). Therefore, a 
similar contribution to the Fermi energy, ep, should be expected. This non-perturbative 
contribution to the Fermi energy induces a modification in the normalization condition of 
the semiclassical density of eigenvalues, eq.(2.13), which should be compensated by the 
normalization of 5u 

X = J = J ^e F -U FP {\)+M. (3.18) 

If ep — ep KB oc /r 6 / 5 exp(-S fec ), then 5 (f(d\/27i)^/e F - U FP ) oc (3~ l exp(— She), and the 
normalization of 5u behaves as expected. Therefore, we obtain that 8u oc exp(— Si, c )5(\ — 
q m ), and the dominant non-perturbative effects in the matrix model free-energy, within the 
quantum mechanical formalism, are 

5F MM oc (3 2 Af J dq >/U F p - e F oc S bc e- Sb % (3.19) 
which are finite, universal, and consistent with the size of the ambiguities of pure gravity. 

4. DISCUSSION AND CONCLUSSIONS 

As we said in the introduction, we should expect to have many different models with 
the same genus expansion as pure gravity or any other fc-even multicritical model, and it 
is necessary to have a criterion to decide which one is describing 2D quantum gravity at 
the non-perturbative level. It seems very difficult to obtain such a requirement from the 
matrix model, where all the results are linked to the double-scaling limit. In fact, even the 
KdV structure, which is a very strong and important result, holds only genus by genus, and 
the matrix models themselves provide examples where it is broken at the non-perturbative 
level. In particular, we shall briefly discuss the solutions of the Hermitian and anti-Hermitian 
one-matrix models with generic potentials, where this is the case. 

The solution of the Hermitian one-matrix models with generic potential, in the one-arc 
sector, is well-known 31 . There exist two functions, — /# ± gg, that satisfy the string 

/ This last factor is apparently missing from ref.[18]. 



14 



equation, eq.(l.l), and in terms of which the KdV flows can be constructed. Their sum, 
fa, is the specific heat of the matrix model, and gn is a non-perturbative function whose 
normalization should be related to the couplings of the odd powers of the matrix in the 
potential. Therefore, the KdV flow structure is not realized in terms of the specific heat 
because of non-perturbative contributions. In the particular case of pure gravity, k = 2, the 
two functions satisfy the Painleve I equation, eq.(1.3), which can be written in terms 
of the specific heat as 

= 2f H9H - \g" H 

3 (4.1) 

The dominant behaviour for large positive values of T is fjj ~ VT, as required, and gu oc 
T~» exp(— iv^y 5 / 4 ). Therefore, the specific heat satisfies an equation different from the 
Painleve I equation, and the difference is non-perturbative. Of course, this provides an 
additional definition of pure gravity to those described in the introduction, but we are again 
restricted to the real solutions of the Painleve I equation, and the problems pointed out by 
David in ref. [7] remain. 

Another interesting example is the solution of the ant i- Hermit ian models 32 . In this 
case, also in the one-arc sector with a generic potential, there is one complex function, 
XA = J A + igA satisfying the string equation, and in terms of which the KdV flows can be 
expressed. The specific heat is given by the real part of this function, and it is obviously 
real, while the imaginary part is non-perturbative. The correlation functions are given by 
the real part of the KdV flow equations 

^ ee (PPO k ) = Re (R' k ) . (4.2) 

Therefore, again, the KdV flow structure is not realized in terms of the specific heat because 
of non-perturbative contributions. In particular, the equations defining pure gravity are now 

= 2f A g A - \g" A 

i ! (4 - 3) 

T = Ia ~ 3 Ia ~ 9a- 

The definition of pure gravity from the anti-Hermitian matrix model has some nice 
relevant features. It was shown in ref. [7] that no real solution of the Painleve I equation can 
be consistent with the loop equations of the matrix model, because all of them have poles 
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along the real axis. Nevertheless, this is not the case with the complex solutions, and there 
is a unique one without poles along the real axis and consistent with the required asymptotic 
behaviour of the specific heat, Ja ~ VT when T — > +00: the "triply truncated solution" 7 ' 10 . 
In contrast with the case of Hermitian matrix models, the complex solutions of the Painleve 
equation are natural in the context of the anti-Hermitian models. Therefore, this definition 
of pure gravity using the "triply truncated solution" is expected to be consistent with the 
matrix model loop equations, and with the reality of the specific heat. It is also worth 
noticing the similarity of this solution to the one proposed in ref.[19] (compare the "Re f" 
in fig.(l) of ref. [11] with fig.(l) of ref.[19]). In fact, their asymptotic behaviours are the 
same also for the large negative values of T, where j a ~ 0, and it could be interesting to 
investigate their relationship further. 

To summarize, we have shown that the stochastic, quantum, or d — 1-like definition of the 
d — matrix models preserves universality in the double-scaling limit. Moreover, the domi- 
nant non-perturbative contributions of the even-A; models are expected to be related to the 
existence of bounces, showing their instability. Nevertheless, we would like to finish by just 
repeating that the essential problem, from the 2D gravity point of view, is how to choose be- 
tween all the possible definitions compatible with the genus expansion of the matrix models, 
and that some new information about the non-perturbative behaviour is needed to solve it. 
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